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ABSTRACT 

A general superconvergence result of conforming finite element method for the elliptic problem is established by 

using L^-projection method. Regularity assumption for the elliptic problem with regular partitions are required. Numerical 
experiment are given to verify the theoretical results. 
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1. INTRODUCTION 

The finite element method (FEM) is a numerical technique for solving problems which are described by partial 
differential equation or can be formulated as functional minimization. To summarize in general terms how FEM works we 
list main steps of the finite element solution procedure below. The first step is to divide a solution region in to finite 
elements. The finite element mesh is typically generated by a preprocess or program. The description of mesh consists of 
several arrays main of which are nodal coordinates and element connectivities. The second step selection interpolation 
function. The interpolation functions are used to interpolate the field variables over the element. Often polynomial depends 
on the number ofnodes assigned to the element. The third step is to find the element properties. The matrix equation for the 
finite element should be established which relates the nodal values of the unknown function to other parameters. The fourth 
step to find the global equation for the whole solution region we must assemble all the element equations. In other words 
we must combine local element equations for all elements used for discretization. Element connectivities are used for the 
assembly process. Before solution 

Boundary condition (which are not accounted in element equations) should be imposed. The fifth step is to solve 
the global equation system. The finite element global equation system is typically sparse, symmetric and positive definite. 
The sixth step is to compute additional results. In many cases we need to calculate additional parameters, for example, in 
mechanical problems strains and stresses are of interest in addition to displacements, which are obtained after solution of 
the global equation system. 

The method was first developed in 1956 for the analysis of aircraft structural problems. Thereafter, within a 
decade, the potentialities of the method for the solution of different types of applied science and engineering problems 
were recognized over the years, the finite element technique has been so well established that today it is considered to be 
one of the best method for solving a wide variety of practical problems efficiently. In fact, the method has become one of 
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the active research areas for applied mathematicians. One of the main reasons for the popularity of the method in different 
fields of engineering is that once a general computer program is written it can be used for the solution of any problem 
simply by changing the input data. Now we wall talk about our main subject. The superconvergence of finite element 
solutions is an interesting and useful phenomenon in scientific computing of real world problems. In recent years, 
supperconvergence for finite element solutions has been an active research area in numerical analysis. The main objective 
in the supperconvergence study is to improve the existing approximation accuracy by applying certain post processing 

techniques that are easy to implement. The superconvergence is obtained by applying the L^-projection method for the 
finite element approximations and their close relatives. In this article, we consider the poisson problem and shall develop a 
general superconvergence results for it is conforming finite element approximation. 

The paper is organized as follows. In section 2, we present review for the elliptic problem and the method to find 
a linear system. In section 3, we develop a general theory of superconvergene result by following the idea presented in 
Wang [3]. In section 4, several numerical experiments are present to support the theoretical result. 



2. PRELIMINARIES AND NOTATION 



Let us consider the poisson problem with Robin boundary condition. Let Q. be bounded and open subset of R 
with Lipschitz continuous boundary (?fl. 



The poisson problem seeks an known function 'If satisfying 



—Au = fin n 



(2.1) 



— hu = gmdil 



(2.2) 



Where T and ft are given function. Here n is the out word normal vector and" - — . 71 where is the 
gradient of u and _i is the Laplacian operator. 



Let 




space of square integrable function define on £2 with inner product and norm 



l 




Further, let be HillbertSobolev space of order 7¥l, that is 



ff™(Q) ={vEL 2 (ft) 



D a V G L 2 (Tt), \a\ < 771 }With inner product 




i 




Iff Ism Q. 
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The weak formulation of (2.1)-(2.2) is to seekUEH ^(Tl)such that 

a(u, v)= (f, i?) n + {g, vy aa VveH\tt) (2.3) 

Where 

a (u, v)= (Vu, Vv) a + {u, v] m = f a Vu. Vvdx + J 3a u. v ds (2.4) 

Let Tfj be a finite element partition of the domain £2 with characteristic mesh size h and lp r be the set of all 
the polynomials defined on k, with degree less or equal to r. Define 

V k = { veH 1 {ti) f v\k£P 1 {k} l Vk€T h } 

The finite element approximation problem for (2.1)-(2.2) is to find'lijj £ Vfy such that 

a (u hf x?)= (f, v) a + is, i?) 9Q yveV h (2.5) 

Where 

a(u hl v) = (Vu hl Vv) a + (wfti = 4 Vu*. Fx? + / aQ u h . v ds (2.6) 

Note that Vfj satisfies the following approximation property. 

^J^-^ll + /l ^-^lll>^ W^Mh^M, ^ r+1 (^) (2-7) 
Define a norm || - || ]_ for 1/ — /f ^ (T?) as follows: 

1MB = a(>, v) = Gfa, Vx?) G + (y, v) da = f a Vv. Vv dx -f- / flQ x? 2 rfs (2.8) 

Definethe basis function (p± £ l 1 ^ with <J?j = Iwhere i=j and (ff j ^ — 0 other wise. 

Our approximate solution lift can be written in terms of it is expansion coefficients and the basis function as 
n 

i=1 

Expand in basis lip, insert into (2.5) and set v=<Pi , i = 1, ... , 71 
n n 

^ c, (T^, V^) n + ^ c ( C^ii P/)«ki = if, Vdo. + {9, Vt)oa 

3=1 i=1 

This is a liner system of equations AX=b with A= (fly} ,X= fCj} ,b= {b^},foi i ,j=l,...,n 
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= (f,<Pi)a + (S,<Pt)an 

Example 2.1: Let the domain O = [0,1] X [0,1] . Also let the exact solution is assumed to be 

w = smtjix) cos(jiy) 

—Ait = 2tt 2 sinOrr) cos (Tiy) in ft 
du 

— + u = g in 3D. 
on 

To find the linear system Ax = b see Figure 1 



n3 < 




> 14 



n2 

Figure 1 

fly = CV^^)fl+ t<Pi,9jXi 

For fej: 

^ 3 = 1 - y, then V<p 3 = 

For &2 ■ 
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then V<p 4 = [_^] 
then SJ<p 2 = [ 1] 



<p 4 = * ~ y, 
<p z =l-x, 
<p 2 =y, thenV<p 2 = [^] 
(Vqt lr V<pO = Vft. 7<Pi dfe x = ^ . . = 
(V<p 2f V<p 2 ) = / fcl V<p 2 .V<p 2 dk x + f k ^ V<p 2 . V<p 2 dfc 2 

i rii rii . 1 roi rQ' 



-iH-El-iGMSl- 



Then 



^1 = 
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= 5lil-xfdx+£y 2 dy = \ 

l 



Then 
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A,= 



■ 2 11 

366 ° 

12 1 

— 0 - 
63 6 

1 21 

6 ° 36 

112 

0 

663 



Since A = A ± + A 2 
Then 



r s-i-i 

o 

3 3 3 

— 15 -1 
Cl 

3 3 3 

-1 S-l 

— 0 

3 3 3 

-1 -15 

0 

3 3 3 J 



A = 



Now to find the load vector 

„ \n cos(nx) cos (ny) 1 
Vu= . \ \ . ' A 
l—nsm{jix) sin Or}') J 

Then / = 2n 2 sin(nx)cos (ny) 



i i 



b{ = (J, <pi) = j /ipidfci = j j In 2 sin (tot) cos(hj) (y — x)dy dx 



0 x 



_4 

7i 



b; = 



f f -7T -8+TT 2 

= (f,<Pi) = j f92dk 1 + J f<p 2 dk 2 = — + - 



2n 
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-Projection Methods 



b* = 



71 

—K -S + 7T' 
2 27F 

2tt 2 
4 



Vit. n 2 + u 



Vit. n 3 H- u 



(— tc sin (rcy) + cos (Try) sin (Try) 
sin(?rr) (irsm(rry) + cosijzy)) 
cos (ttj) (sin(Trx) — ncos(jzx)) 
cos (Try) (ttcos(tt?c) + sin(nx)) 



1 1 

= j -sin(nx) (1 — x) (te + J —Ticos (7ry}t£y 



Then J&* = 



1 
1 
1 



= fcj theji x = A 1 b 

Lemma 2.1 

There exists a constant, C, independent of tl such th 

\a(v f w)\ < CIMlillwlli, Vw^H 1 ^) 

Proof: since the problem (2.1) - (2.2) satisfy the weak form 
Where 
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1 Vw. Vi? dx 


+ 


1 w.vds 









a(w,v} = (Vw, Vv) a + (w, i?) 3 q = I Vw. Vt? dx + I iv. i? 

By using Schwarz inequality we obtain 

|a(w,t?)| = 

^ l|Vw|| i -(rolhMlj?"(H> + llwll^a^lMl^aQ) 

By the trace inequality there is a constant, C, such that 

IMIl 2 oq) ^ c IML 

We have 

|a(w,v)| < IIVwIU^IIVvH^gd + C HwlliClblli 

Since 1 1 VlV || L z (Q> < llwllj 

l l 
WhereUVwIL^ = (f fl \Vw\ 2 dxf ,\M\l= (f n \Vw\ 2 dx + f^ds)* 

Then we have 

k(w,iOI< DwDJHIi + ^DwDilklli 

<L0-+C*)\M\i\\v\\i 
< CllwII Jblli 

Where C = 1 + C 2 

And thus the proof is complete. H 
Theorem 2.1 

Let tt and lift, be the solution of (2.1)- (2. 5) respectively. Then, there exist a constant, C f independent of 
h such that 

\\u-u h \\ 1 ^Ch k \\v\\ k+1 _ 

Proof: 

Since 
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Hu-UftHi < inf ||u- f/lli 

and by using (2.7) then 

Hu-uJa < Ch k+1 \\u\\ k+1 forU€H^Hm 

3. SUPERCONVERGENCE BY £ 2 -PROJECTION 

-projection technique was introduced by Wang [3]. It projects the approximate solution lift to another 
finite element dimensional space associated with a coarse mesh. The difference in size of the two meshes can be used to 

achieve a supperconvergence. Now, we start with defining a coarse mesh 7^ where T » h satisfying: 

T = k a (3.1) 

With CL £ £0,1}. For now, the parameter Ct will play an important role in the post processing. Define a finite 
element space 

Herein consists of piecewise polynomials of degree T associated with the partition J^-. 

Let Q?be the L^-projector onto finite element space V^-. Ql^Can be considered a linear operator (projection) 
from I? (ft) onto the finite element space l^.We will assume an H s , S ^ 1, regularity for the solution of 

(2.1), (2.2): 

II « ll s < C || / |[ s _ 2 

(3.2) 

Lemma 3.1 

Suppose that (3.2) holds with 1 S 5= It -+- 1 and the finite element space 1^ C fi s ^ Then, there 
exists a constant, C 7 independent of h and T, such that 

|| Q^U - Q T U h || < Ch'-^^^WU-U^, (3.3) 
Where ^^(O.l^s the parameter defined in (3.1). 

Proof: The definition of ||. ||and Q. z give 

|| Q T u - Q T u h ||= SUp || (Q^U- Q T U ht <p)|| 

From definition of the L^-projection for Q T , we have 

(Q T u - Q. T u hf <p) = (u- u hf Q?<pl 
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Then 

II Q^u-Q T u k ||= snp vefi3( flj (l(pl=1 ||Cu-« fcj Q T <p)\\ , (3.4) 

Consider the following problem: find W€fi * £flj such that 

a (w, v) = (Q T <p, vyVvcH 1 (fi) 0.5) 

It's follows from (3.2),wff s (fi) f\ H 1 (fl)and 

||W|| S < lle^lL-z, (3.6) 
For some constant C Using the inverse inequality on the right hand side of (3.6), we have. 

\\w\\ s <CT^ 2 -^\\<p\\o (3.7) 
Subtracting (2.5) from (2.3) gives 

a{u- U kl v) = 0, VV€V h (3.8) 
By replacing V by It It in (3.5) and using (3.8), we obtain 

(u - u hl Q T <p) = a(u- u hl w) 

= a{u - U hl W — 1?) (3.9) 
Where V€V h .Using Lemma (2.1) and Theorem (2.1) we obtain from (3.9) 

\(u-u h ,Q T <p)\ = a\{u-u h ,w- v)\ 

< C inf llu-UftllJIw-T?!!! 

^Ch^Wu-uMMU. am) 

From (3.6), (3.7) and (3.1) we obtain 

\(u - u h , Q T <p)\ < Ch?' 1 || Q T <p || s _ 2 || u - u h || x 
< ck^H" 11 ^ 0 ' 2 -^ \\u-u h || -l || <p || 

Then 

|| Gr^ft II < C^- 1+minfe2 - s} || U- IU.. 
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Theorem 3.1 

Assume that (3.2) with 1 < s < k + 1 and the finite element space V^. C ^ (Tl}lf the exact solution 
U G H k+1 (fi) n H r+1 (Ti) n H 1 (fl), then there exists a constant, C f such that 

||u - Q x u h \\ 4- fc*|[Vr(u - Q xUh )\\ 

< Ck« ( r+ti\\ u \\ r+1 +C^|t„-« ll || 1 (3 .10) 

Where £7 = S ~ 1 H" Ct 111111 (0,2 — S) and U h is the finite element approxima-tion of the solutionis. 
Proof: Since Q^is the L^-projection, then 

\\U- Q T U\\ * CT r+1 ||u|| rHhl < Cft ff(r+l3 ||u|| r+1 . (3.11) 
Since 

||it - Q T « A || < ||« - Q T u|| + ||Q T « - Q x u h \\ 

Combining (3.3) and (3.11) in above inequality gives 

\\u-Q r u h \\ < Ch* ( r + ti\\u\\ r+1 +Ck'- x+aa **& 2 -*)\\u-vift if 

Which complete the estimate for \\tL Q T Ufa || in (3.10) 

Now we estimate || V T (VL — Q^lifr || as follows. Since, Q T is L^-projection and using (3.1) we have 
||V T (M " Q x U h )\\ < CT r ||M|| r+ i ^ Cft"-||tt|| r+1 . (3.12) 
Multiplying both sides of (3.12) by fl X , we get 

h*\\V T (u-Q T u h \\ < £7i*< r+1 )||ii|| r+1 (3.13) 

The inverse inequality gives: 

IIV r «? T u - Q 7 u h )\\ < Ct^ICM* - Q T u h \\ (3 i4) 

Using (3.1) and multiplying (3.14) by h a the above inequality implies 

h a \\V x (Q v u- Q T u h )\\ < C\\Q x u-Q T u h \\ 

< ch s - 1+<xmlR < a2 -*> \\u - uJk, (3.15) 

and 

h*\\V T {u - Q T uJ|| < fc"||V T (u- Q t m)|| + A«|[V, CQ T «- Q T « h )ll 
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Then 

II (« - Q x uO II II V T (u- Q T UjJ || 

< Cfo ff(r+1] || u || r+1 + C?i fc+s - 1+amin ^ 2 -^|| u || / . + 1 

The above error estimate is optimized if £t is selected so that 

a(r + 1) = k+ s — 1 + a min(0,2 - s). 

Solving CC from above yields 
kA-s-1 

CL = ; — ; (3.16) 

r + I— rain (0,2— s) 

The corresponding error estimate is given by 

II (u - Q x u£ || +h a || V T (u- Q T u h ) || 

<Ch^ +1 \\\u\\ r+1 +\\u\\ k+1 ) (3.17) 

Theorem 3.2 

Suppose that (3.2) holds with 1 $ + 1 and the finite element space be the 

finite element approximation of the solution. Then the post-processed of U is estimated by (3.17) with 

k+s-1 



a = 



r + 1 - min(0,2 — s) 

4. NUMERICAL EXAMPLES FOR CONFORMING FEM BY £ f 2 -PROJECTION METHOD 

In this section, we present several numerical examples to verify Theorem3.3. The triangulation of 7\ is 
constructed by: (1) dividing the domain into an 71 X 71 rectangular mesh and (2) connecting the diagonal line with the 

positive slop. Denote h, = — as the mesh size. The finite element spaces are defined by 

■n 3 

V T = {v€L 2 m:v\keP z (K),Vk6T T } 

Example 4.1: Let the domain fi = [0,1] X [0,1] and the exact solution is assumed as 

u = (r 3 — x 2 ') (y 3 — y 2 ') 

Table 1 shows that after post-processing method, all the errors are reduced. The exact solution in L^-norm of 
\\u — Qx^h Hi has the similar convergence 
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Rate as ||u — ||. There is no improvement for the U in L^-norm. However, the error in the M ^-norm has a 

higher convergence rate, which is shown as Q^tl^"^ for ||V T (u — Q T Zt h ) ||, see Figure 2(a) and (b). Figure 3(a) and 
(b) give the results for the finite element approximation of (2.1) and (2.2) before and after post-processing. 



Table 1: Error on Uniform Triangular Meshes 7}, And T T 



h 


| U — Uk | i 


. it - Uv | |K* — Q T U h | -t | U - Q T U h \_ 


V 


0.8936e-2 


0.6620e-3 


0.6455e-2 


0.6340e-3 


3" 3 


0.2842e-2 


0.5986e-4 


0.1172e-2 


0.5707e-4 


4" 3 


0.1212e-2 


0.1066e-4 


0.3510e-3 


0.1015e-4 


5" 3 


0.6224e-3 


0.2795e-5 


0.1404e-3 


0.2661e-5 


6" 3 


0.3604e-3 


0.9358e-6 


0.6698e-4 


0.8909e-6 


o(k k lk= 


0.9752 


1.9915 


1.3873 


1.9933 



10" 
10" : 
10- 

itf 

10"' 



-a-l|U-U h l| 
-+-I|u-Qu h || 

»Q(h 2 ) 



10" 



10" 



10 10 

(a) 



10° 



10"' 









-h-Iu-Qu^ 


v 


^O(h) 









10" 



10" 



10' 



10° 



(fe) 



(a) (b) 
Figure 2: (a) Convergence Rate of L J -Norm Error; (b) Convergence Rate of ff 1 -Norm Error 




Figure 3: Result for U = (x 3 - X 2 }(y 3 ~ y 2 ) In Example 4.1. (a) Surface Plot of 
Approximation Solution (b)Surface Plot of Approximation Solution Q^t^ 

Example 4.2: Let the domain ft = [0,1] x [0,1] 

and the exact solution is assumed as 

u = x 5 cos (y) 

From the results shown in Table 2, it is clear that the exact solution U in H *-norm has the superconvergence but 
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there is no improvement in L^-norm, see Figures 4 (a) and (b). The finite element solution given in Figures 5 (a) and 
(b).This agree well with the theory. 

Table 2: Error on Uniform Triangular Meshes and T T 
h | u — Uf.. | 1 | it — Ui., | | u — Q T u h | ! | u — Q T u h \\ 

3" 3 0.3407e-10.5738e-3 0.8608e-20.4887e-3 
4" a 0.1441e-10.1023e-3 0.2618e-20.8709e-4 
5~ 3 0.7384e-20.2684e-40.1065e-20.2288e-4 
6~ a 0.4274e-20.8991e-50.5135e-30.7648e-5 
0 (k k ), k= 0.9982 1 .9986 1 .35561 .9993 




(a) (b) 

Figure 4: (a) Convergence Rate of Z/-Norm Error; (b) Convergence Rate of ff*-Norm Error 




Figure 5: Result for u = % 5 co-s (y) in Example 4.1 (a) Surface Plot of Approximation 
Solution (b)Surface Plot of Approximation Solution Q T U h 
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